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The mass spectrum of baryons in the spin-3/2 sector is computed in quenched lattice QCD using 
a tadpole-improved anisotropic action. Both isospin 1/2 and 3/2 (the traditional decuplet) are 
considered, as well as members that contain strange quarks. States with positive and negative 
parities are isolated by parity projection, while states with spin-3/2 and spin-1/2 are separated by 
spin projection. The extent to which spin projection is needed is examined. The issue of optimal 
interpolating field is also investigated. The results are discussed in relation to previous calculations 
and experiment. 

PACS numbers: 12.38.Gc, 14.20.Dh, 14.20.Gk, 14.20.Jn 



I. INTRODUCTION 

The mass spectrum of hadrons represents a fundamen- 
tal manifestation of the long-distance dynamics of quarks 
and gluons as governed by QCD. Non-perturbative cal- 
culations through numerical simulations on a space-time 
lattice provide a method to obtain this spectrum from 
first principles. The computation of the hadron spec- 
trum using lattice QCD started in the early 1980's [llQl- 
The modern era in lattice QCD calculation of the hadron 
spectrum started with the results of the GFll group ^. 
The benchmark calculation of the quenched light hadron 
spectrum using the standard Wilson action has been per- 
formed by the CP-PACS collaboration |3,I3- For state- 
of-the-art computations using dynamical configurations 
that involve the study of baryon mass spectrum, see, for 
example, ,6„ uL iSb iSj ■ For reviews on baryon mass spec- 
trum, see, for example, p!ol.[ll|. 

Most of the lattice computation of the light hadron 
spectrum has been limited to the ground states. It is im- 
portant to extend the successes beyond the ground state. 
The rich structure of the excited baryon spectrum, as 
tabulated by the Particle Data Group 12], provides a fer- 
tile ground for exploring how the internal degrees of free- 
dom in the nucleon are excited and how QCD works in a 
wider context. One example is the parity splitting in the 
low-lying N* spectrum. The nucleon A^(938) has positive 
parity, while its negative parity excitation, S'ii(1535), has 
a much higher mass. The spontaneous chiral symmetry 
breaking in QCD is thought to be responsible for the 
splitting. Without it, QCD would predict exact parity 
doubling in the baryon spectrum. The study of the ex- 
cited mass spectrum is a critical part of the experimental 
program at Jefferson Lab. Lattice QCD has a number of 
advantages in helping understand the N* spectrum. One 
can systematically study the spectrum sector by sector, 
with the ability to dial the quark masses, to separate the 
parities exactly, to project out the spin components, and 
eventually to dissect the degrees of freedom in the QCD 
vacuum most responsible for the spectrum. There have 
been a number of lattice calculations of the N* in the 



spin-1/2 sector [H El El [H IH 111 El H IH IH HI . 

In this work, we focus on the spin-3/2 sector. In ad- 
dition to the usual baryon decuplet with isospin-3/2 and 
spin-3/2, we study the isospin-1/2 and spin-3/2 fam- 
ily which has only received limited attention so far. A 
preliminary study was reported in Ref. |24j . A calcula- 
tion using the FLIC fermion by the Adelaide group was 
done in Ref. (23 ■ Other mehtods for constructing higher 
spin states have been proposed by the LHPC collabora- 
tion [21123,123 ■ Here, we use a different interpolating 
field as the one used in Ref. [23 . We also extend the cal- 
culation to include states that contain the strange quark. 
The goal is to establish the basic features in terms of 
spin-parity on the lattice. 



II. CALCULATION DETAILS 

Excited states composed of light constituents are both 
large in size and mass. Their calculation imposes se- 
vere signal-to-noise problems. The use of an anisotropic 
lattice can help alleviate the problem. A fine lattice in 
the temporal direction enables the correlator to be ob- 
served over many time slices at short separations, while 
the coarse spatial spacing allows large spatial volumes 
the states demand. We use the anisotropic gauge action 
given in Ref. |29| : 
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where u., — (iReTrP, 



^/■* is the spatial tadpole factor, 
Pss' denoting the spatial plaquette. Ut is the temporal 
tadpole factor, we set Uj = 1 in this simulation, flc = 
J2c ^^eTi{l — Wc), with Wc denoting the path-ordered 
product of link variables along a closed contour C on the 
lattice, rigp includes the sum over all spatial plaquettes 
on the lattice, iltp indicates the temporal plaquettes, fisr 
denotes the product of link variables about planar 2x1 
spatial rectangular loops, and figtr refers to the short 
temporal rectangles (one temporal link, two spatial). 



For the quarks, the anisotropic D234 action of Ref. |30|, |31| is used with the following Dirac operator, 
M^,34 = m(l + \ram) + E{7.AW _ ^^^a(3) + ,[_« a(2) _ « ^ .^^p,. + |!a(4)]}. 
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Here AJT^ is the n*'' order lattice covariant derivative, A^^) = A^^^A^^) = A^^^A^^), and A^") = A(2)a(2). The terms 
proportional to r are generated by a field redefinition and thus represent a redundant operator. With the help of the 
gauge-covariant first- and second-order lattice derivatives. 
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and 
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the action can be cast into the standard form of 
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With the specific choice of the factors b^ = ^, Cf^ = 
^^, and r = 2/3, and an improved version of the field- 



strength operator free of 0{a ) errors, F^u 
relative 0{a^) correction. 



(clover) 
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the action has only 0{af,a'l) classical errors. This ac- 
tion consists of three types of interaction terms: nearest- 
neighbor, next-nearest-neighbor, and the clover term. 
Both gauge action and quark action have tadpole- 
improved tree-level coefficients to reduce unwanted quan- 
tum fluctuations. The hopping parameter k is related to 
the bare parameters by 
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In this calculation, we use an 10'^ x 30 anisotropic lat- 
tice with anisotropy ^ = as /at = 3. The spatial lat- 



tice spacing a^ w 0.24 fm determined from the Som- 
mer scale rp. The lattice coupling (3 — 2.4. In 
all, 100 configurations are analyzed. On each con- 
figuration 9 quark propagators are computed using a 
multi-mass solver, with quark masses ranging from ap- 
proximately 780 to 90 MeV. The nine k values are: 
Ki_9 = 0.30,0.31,0.32,0.33,0.34,0.345,0.35,0.355,0.36 
They correspond to pion mass in the range of 2.11 to 
0.68 GeV, and the mass ratio -k/ p from 0.95 to 0.65. 
The strange quark mass corresponds to the seventh kappa 
value (k — 0.350). The critical kappa value determined 
from TO^ is Kc = 0.3705(3). The source is located at 
{x,y,z,t) = (1,1,1,2). We use Dirichlet boundary con- 
ditions in the time direction. 

A Gaussian-shaped, gauge-invariant smearing func- 
tion [33 in spatial directions. 



a J2 (^»' 

i=l,2,3 



uHx)) 



N 



(8) 



was applied both at the source and at the sink to in- 
crease the overlap with the states in question. So for a 
given interpolating field operator, one can construct four 
types of correlation functions with the source-sink combi- 
nations of local-local (LL), smear-local (SL), local-smear 
(LS), and smear-smear (SS). In Eq.®, a is the couphng 
strength at which the neighboring links are brought in, 
and N is the number of iteration times. We used a — 0.25 
and iV = 10 in all cases. We found that the SL gives the 
best signal so the results presented in this work are from 
this combination. 

We consider the full interpolating field with the quan- 
tum numbers I{J^) = 5 (| ) ^^ proposed in Ref. [33, 



X^ = eabc (u'^^C^.ipS') (g^-P - \l''l'\ 15U'- 
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It satisfies the condition 7^x^ = for spin-3/2 fields. 
The superscript T denotes transpose. The C = 7472 
is the charge conjugation matrix. The Dirac 7 ma- 
trices are Hermitian and satisfy {71/, 7r} = '^^ut, with 
ai^r = TJ^[7 t/, ^ tI- We follow the gamma- matrix notation 
of Sakurai p34| . The /x, p are Lorentz indices and summa- 
tion over p is implied. The antisymmetric eatc ensures 
that the state is color-singlet. The interpolating fields of 
the spin-3/2 E* and S* are obtained by properly chang- 
ing the quark field operators. For example, one can get 
the spin-3/2 S* state interpolating field by substituting d 
with s; and S* by replacing u with s. This interpolating 
field has 5 terms as compared to the standard interpo- 
lating field for the nucleon, 



(10) 



X{x) = <^abc u"- {x)C^5d (x) u''{x) 



So computationally, it is 25 times more expensive. Fur- 
thermore, since the full 4x4 matrix in Dirac space (as op- 
posed to only diagonal elements) is needed to carry out 
the spin projection described below, an extra factor of 4 
is needed, making this interpolating field 100 times more 
expensive than a standard nucleon mass calculation. 

Despite having an explicit parity by construction, the 
interpolating field couples to both positive and negative 
parity states. A parity projection is needed to separate 
the two. In the large Euclidean time limit, the two-point 
correlation function with Dirichlct boundary condition 
in the time direction and at zero spatial momentum be- 
comes 
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where f^i, is a function common to both terms. The rel- 
ative sign in front of 74 provides the solution: by taking 
the trace of G^,y(i) with (1 ± 74)/4, one can isolate M+ 
and M_, respectively. 



The interpolating field in Eq. © couples to both spin- 
3/2 and spin-1/2 states. To project a pure spin-3/2 state 
from the correlation function Gp,i,, we use a spin-3/2 pro- 



jection operator |25ll3 
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P^^(3/2) =ffM''--7M7i^-^(7-P7MP^+PA'7>^7-p) (12) 

The corresponding spin-1/2 state can be projected by 
applying the projection operator 



P^,(l/2)=5^,-P^,(3/2) 
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The projection is done after the correlation functions are 
generated, with no need to generate new quark propaga- 
tors at the source. Only zero spatial momentum {p = 0) 
is considered in the projector. To use this operator and 
retain all Lorentz components, one must calculate the 
full 4x4 matrix in Dirac and Lorentz space of Gp„{t). 
Using the projection, we have 
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However, to extract the mass, only one diagonal pair of 

Lorentz indices is needed, reducing the amount of calcu- 

1/2 
lations required by a factor of 4. We calculate G33 and 
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III. ISOSPIN-1/2 AND SPIN-3/2 BARYONS 

Fig-ffldcmonstrates results for the correlation function 
for both parities in the nucleon channel at the smallest 
quark mass considered. In the positive-parity channel, 
one can see that this correlation function is almost com- 
pletely dominated by the 1/2-t- component. Spin projec- 
tion reveals two different exponentials from the spin-3/2 
and spin-1/2 parts, with the spin-3/2 state being heavier 
than the spin-1/2 one (a steeper fall-off), in agreement 
with the ordering in experiment. The expected relation 
in Eq. l(TBjl is indeed satisfied numerically, providing a 
non-trivial check of the calculation. 

A further check of the calculation is provided by the 

1/2 
fact that the mass extracted from G^u {t) is degener- 
ate with that from the conventional G{t) for the nucleon 
ground state using the standard interpolating field in 
Eq. (|10|) . One can see that spin projection is crucial 
in this channel. Without it, one would get a false signal 
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FIG. 1: The various correlation functions (unprojected, spin- 
3/2 projected, spin-1/2 projected) for the nucleon states in 
the positive-parity (top) and negative-parity (bottom) chan- 
nels at the smallest quark mass (k = 0.36). 



for spin-3/2 since it is dominated by the spin-1/2 com- 
ponent state. The large error bars is a sign of sensitive 
cancellations in the projection procedure. 

The situation in the negative-parity channel is oppo- 
site, as shown in Fig. ^ Here the signal is dominated by 
the 3/2— state, so one would get a spin-3/2 signal with- 
out spin projection. The results also show a similar fall- 
off for the 1/2— state and the 3/2— state, in accord with 
the experimental states of iV*(1535)i" and A^*(1520)f " 
which are close to each other. We checked that the con- 
dition in Eq. (|16|l is also satisfied. According to Fig. ^ 
we can get a rough idea about the range of time slices 
we should choose to extract the baryon masses. For ex- 
ample, to get good fitting result of iV*(3/2-|-), we should 
choose time slices earlier than 12, however, the mass of 
A^(l/2-|-) can be extracted from much later time slice. 
To find out more specific fit time window, we use effec- 
tive masses and extract the baryon masses from a plateau 
area. Fig. |21 presents the effective masses in the spin-3/2 
sector at four quark masses that correspond to the heav- 
iest quark mass, and two quark masses in the middle, 
and the lightest quark mass). For iV*(3/2+) the sig- 
nal is weak because the correlation function for positive- 
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FIG. 2: Effective mass plot for the iV*(3/2-|-) state (top) and 
N*{3/2—) state (bottom) at selected quark masses. 



parity is dominated by the N*{l/2+) state. Only a rough 
plateau from time slice 7 to 9 can be found and the mass 
of iV*(3/2-|-) is extracted from this time window so the 
results for this state should be taken with caution. To 
access later time slices, a large number of configurations 
are needed to increase the signal-to-noise ratio. The sig- 
nal for 7V*(3/2— ), on the other hand, is much stronger 
since it is the dominant component in the negative par- 
ity channel. A nice fiat area can be found between time 
slices 8 and 12. The mass of the iV*(3/2— ) is extracted 
from time slice 10 to 12. 

Fig. O displays the effective masses in the spin-1/2 sec- 
tor. Here A^(l/2-|-) is the dominant component, while 
A'^(l/2— ) is the weaker one, in their respective parity 
channels. The N{l/2+) state is extracted from time slice 
11 - 14. For A^(l/2— ), a rough flat area can be found 
from time slice 7 to 12, but we use 9 to 11 since the value 
from later time slice has smaller systematic error. To 
show the quality of the nucleon from the two different 
measurements, a comparison of the effective masses of 
the nucleon at two kappa values is given in Fig. ^ The 
projected nucleon at the smallest pion mass shows some 
instability beyond time slice 15. 

Fig. |31 presents results of the mass ratios extracted 
from the correlation functions for the N* states to the 
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FIG. 3: Effective mass plot for the A''* (1/2+) state (top) and 
N*{l/2~) state (bottom) at selected quark masses. 
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FIG. 4: Comparison of nucleon effective masses at the small- 
est and largest pion masses from two ways: one from the 
standard operator in Ea.f ll()|l . the other from the 1/2+ com- 
ponent of the full operator in Eq.(|5J. 



nucleon ground state as a function of the mass ratio 
(tt/p)^. We use mass ratios because they have minimal 
dependence on the uncertainties in determining the scale 
and the quark masses, and have smaller statistical er- 
rors. In this figure, we slightly shift the points for the 
N*{l/2—) state to the right hand side to avoid overlap- 
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FIG. 5: Mass ratio of the projected A'^* states to the ground- 
state nucleon as a function of the mass ratio squared (j^ j pf' . 
The four lattice states with separated spin-parity are symbol- 
coded (and color-coded if color is visible) with the experimen- 
tal candidates which are indicated on the left at the physical 
point. 



ping. The mass of the nucleon ground state we use in 
the mass ratios is produced on the same configurations 
from the standard interpolating field. It is encouraging 
to see that the 1/2+ nucleon obtained from the standard 
operator agrees with that from the projected 1/2+ nu- 
cleon in the spin-3/2 operator (as indicated by the ratio 
of 1). Clear splitting is seen between the 7V(l/2+) state 
and its parity partner A^*(l/2— ) state in this figure. The 
two negative-parity states iV*(l/2— ) and A^*(3/2— ) are 
degenerate within errors, which is consistent with the 
fact that the corresponding states A^(1535)l/2— ) and 
A^(1520)3/2— in the observed spectrum are very close 
to each other. We do similar calculations in the E* and 
S* channels. The results for the S* channel is shown in 
Fig. In the PDG jljl the two positive-parity states 
i;(1193)(l/2+) and I](1385)(3/2+) are weh established, 
which we identify as S(l/2+) and I]*(3/2+) on the lat- 
tice, respectively. The situation in the negative-parity is 
unclear. The two possible candidates for the lattice states 
i;*(l/2-) and I]*(3/2-) are 2-star states i;(1620)l/2- 
and I](1580)3/2— , respectively. In the figure, they have 
question marks, and the four lattice states are matched 
one to one with the experimental candidates. 



A similar situation exists in the a^ channel shown in 
Fig. El The two well-established positive-parity states 
S(1318)l/2+ and S(1530)3/2+ are identified with their 
lattice counterparts. The two possible candidates for the 
negative-parity lattice states S*(l/2— ) and S*(3/2— ) are 
3-star S(1690) state with unknown spin-parity and the 
3-star S(1820)3/2— state, respectively. All of the mass 
ratios in this sector are listed in Table Ulll 



TABLE I: Mass ratios for the four states in the isospin-1/2 and spin-3/2 family after parity and spin projections. The results 
are from the full interpolating field in Eq. ^ . The numbers in brackets are statistical errors in the last digits. The last column 
indicates the time window from which the results are extracted. The absolute values for the pion mass and the nucleon mass 
are also given for conversion purposes. 



m^ (GeV) 


2.10(1) 


1.87(1) 


1.65(1) 


1.43(1) 


1.20(1) 


1.08(1) 


0.96(1) 


0.82(1) 


0.68(1) 


10-15 


ruN (GeV) 


3.45(1) 


3.12(1) 


2.80(1) 


2.49(1) 


2.19(1) 


2.04(1) 


1.88(1) 


1.72(2) 


1.56(2) 


10-15 


rrin/mp 


0.948(1) 


0.935(1) 


0.916(2) 


0.889(2) 


0.847(2) 


0.817(3) 


0.777(4) 


0.723(4) 


0.643(5) 


10-15 


N{l/2+)/N 


1.00(0) 


1.00(0) 


1.00(0) 


1.00(0) 


1.00(0) 


1.00(0) 


1.00(0) 


1.00(1) 


1.00(2) 


11-14 


N''{l/2-)/N 


1.14(3) 


1.16(3) 


1.19(4) 


1.22(5) 


1.28(6) 


1.32(8) 


1.36(10) 


1.42(13) 


1.50(17) 


9-11 


N*{i/2+)/N 


1.28(4) 


1.31(4) 


1.36(4) 


1.41(4) 


1.48(5) 


1.52(6) 


1.56(6) 


1.61(7) 


1.67(8) 
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N*{3/2-)/N 


1.15(1) 


1.17(1) 


1.20(1) 


1.24(2) 


1.30(2) 


1.33(2) 


1.38(3) 


1.44(3) 


1.53(4) 


10-12 


E*(l/2-h)/iV 


0.85(0) 


0.87(0) 


0.90(0) 


0.93(0) 


0.97(0) 


1.00(0) 


1.04(1) 


1.09(2) 


1.16(3) 


15-17 


E*(l/2-)/iV 


1.03(6) 


1.07(7) 


1.12(8) 


1.18(10) 


1.27(12) 


1.33(14) 


1.40(15) 


1.49(18) 


1.60(22) 
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0.96(9) 


1.00(9) 


1.05(11) 


1.10(12) 


1.18(14) 


1.22(14) 


1.27(16) 


1.32(17) 


1.38(18) 


10-12 


E*(3/2-)/iV 


1.03(1) 


1.07(1) 


1.12(1) 


1.18(2) 


1.26(2) 


1.32(2) 


1.38(3) 


1.46(3) 


1.56(4) 
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H*(l/2+)/iV 


0.71(0) 


0.75(0) 


0.79(1) 


0.84(1) 


0.91(1) 


0.95(1) 


1.00(1) 


1.05(1) 


1.13(2) 
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H-(l/2-)/iV 


0.89(2) 


0.95(3) 


1.02(3) 


1.11(4) 


1.22(4) 


1.29(5) 


1.38(6) 


1.48(7) 


1.62(8) 
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0.84(7) 


0.88(7) 
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1.07(9) 


1.12(10) 


1.18(10) 


1.26(11) 


1.36(13) 
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1.41(10) 


1.52(12) 


1.66(14) 


13-15 



2.0 

1.8 

. 1.6 



1.4 



S 1 



.2 
1.0 
0.8 



♦ 1620-(?)/93a 
Ol580-(?)/938 

n 1385V938 

• 1193V938 



n 3/2+ 


O 3/2- 


• 1/2 + 


♦ 1/2- 




0.0 0.2 0.4 0.6 

(m /m y 



1.0 



1.4 
1.3 
1.2 



+ 

C\2 1.1 



1.0 
0.9 
0.8 



First 1 



938+7938 



0.0 



0.2 



i 11 A i »4 



0.4 0.6 



0,8 



1.0 



FIG. 6: Same as in Fig. |3 but for the E states. 
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FIG. 7: Same as in Fig. |3 but for the H states. 



FIG. 8: The mass difference between the result from the full 
interpolating field and that from the first term of the inter- 
polating field for the N{l/2+) state. 



IV. THE ISSUE OF INTERPOLATING FIELDS 

The computational cost of evaluating each Lorentz 
combination in Eq. @ is relatively high, about 100 times 
that for the ground-state nucleon. Therefore, it took a 
long time for us to get the 100 configurations. In Ref. |25|. 
only the first term of the full interpolating field, 

Xm = ^abc (w^^CtsTmC?') 75U^ (17) 

was considered. In our calculation, we deliberately sepa- 
rated the contributions from individual terms, so we are 
in a position to investigate possible differences between 
the two interpolating fields. 

We found that the projected spin-3/2 components are 
the same within statistical errors. In the spin-1/2 sec- 
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FIG. 9: The various correlation functions (unprojected, spin- 
3/2 projected, spin-1/2 projected) for the Delta states in the 
positive-parity (top) and negative-parity (bottom) channels 
at the smallest quark mass. 



tor, however, the situation is different. Fig. |S1 shows the 
comparison in the 1/2+ channeL There is a significant 
difference (up to 30%) between these two results. Similar 
difference is found in the 1/2— channel. This result is dif- 
ferent from that in Ref. z5]. So as far as spin-3/2 states 
are concerned, the first term of the interpolating field can 
be used to do the calculation. However, if the masses of 
spin-1/2 states are desired from such interpolating fields 
via spin projection, the full interpolating field is required. 



V. ISOSPIN-3/2 AND SPIN-3/2 BARYONS 

For a A+ state, the interpolating field is 

A+ 1 



Xfi 



V3 



iabMu^''Crd''X- + iu^''Cj^u'')d% (18) 



Interpolating fields for other dccuplet baryons can simi- 
larly be obtained by appropriate substitutions of quark 
fields. 

First of all, we would like to show the correlation func- 
tions since these are the bases of mass extraction. Fig. 
presents the correlation functions for A states in both 
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FIG. 10: Mass ratio of the projected A states to the ground- 
state nucleon as a function of the mass ratio squared (iv/p)'^. 
The four lattice states with separated spin-parity are symbol- 
coded with the experimental candidates which are indicated 
on the left at the physical point. 



parity channels. Again, the relation in Eq. (|16|l is satis- 
fied. The spin-3/2 component almost completely domi- 
nates in both channels. The 1/2+ component is weaker 
by several orders of magnitude. This indicates that the 
interpolating field in Eq. 1181 has a small overlap with 
spin-1/2 states. Nonetheless, wc could extract a discern- 
able signal with 300 configurations. The mass of the spin- 
3/2 states can be extracted from relatively large time 
slices, while the mass of the spin-1/2 states are extracted 
before time slice 10 since the data becomes very noisy be- 
yond that point. Another observation is that A(l/2+) is 
heavier than A(3/2+), while A(l/2— ) is about the same 
as A(3/2-). 

As usual, the fitting is done on the effective masses, 
using the same spin and parity projection techniques. 
Here we show the results directly without showing the 
effective mass plots. Fig. ^] shows the results in the 
A channel. The trend of the A(3/2+)/A^ data points 
with decreasing quark masses is clearly toward the ob- 
served ratio of A(1232)/Af(938). The splitting between 
A(3/2— ) and A(3/2+) is consistent with that in exper- 
iment (470 MeV). There is a small hint that A(3/2-) 
lies above A(l/2— ), but the errors are too big to resolve 
the two states which are close to each other in exper- 
iment. For better viewing, the points for the A(l/2— ) 
state have been slightly shifted to the right hand side, and 
the points for the A(l/2+) state have been shifted to the 
left hand side. The large errors show that this state is 
difficult to extract with limited statistics. The splitting 
between A(l/2+) and its parity partner A(l/2— ) also 
appears consistent with that in experiment (290 MeV). 
The signal of this A(l/2+) state is the weakest in the 
four A states. In all, the splitting pattern of these states 
is consistent with that observed in experiment and with 
.Ref. [23, despite weak signals for the A(l/2±) states. 

Similar patterns exhibit in the E* and S* states in the 
decuplet, as shown in Fig.lllland Fig.^] Strong signals 
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FIG. 11: Same as in Fig. ^1 but for the E states. 
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FIG. 12: Same as in Fig. ^| but for the ^ states. 



and stable results for the spin-3/2 states, relatively weak 
signals and big error bars for the spin-1/2 states. Note 
that there is a systematic improvement in the signal for 
the 1/2+ states as the nmxiber of the strange quarks in- 
creases from in A to 1 in E* to 2 in S*. This can be 
attributed to the stabilizing effects of the heavier strange 
quark. The experimental situation for the S* states is 
not clearly settled yet ^^. That is why we put ques- 
tion marks on some of the states. The 5(1820)(3/2— ) 
state has 3-star status which we identify with our 3/2— 
state on the lattice. There is a 3-star state S(1690) with 
unknown spin-parity. We identify this state with our 
1/2— state on the lattice. Our highest state S*(l/2-|-) is 
identified with the 3-star state of S(1950) which also has 
unknown spin-parity in the Particle Data Group. All of 
our results in this sector are summarized in Table Ivl 



used. We analyzed 100 configurations despite a big in- 
crease in computing demand as compared to a truncated 
version of the interpolating field. Four states with defi- 
nite spin-parity are isolated for each particle type using 
parity projection and spin projection. The need for spin 
projection is clearly demonstrated in the positive-parity 
channel whose correlation function is dominated by the 
spin-1/2 component. Clear signals are obtained for both 
the spin-projected iV*(3/2±) and the N*{l/2±) states, 
although the latter are usually weaker, resulting in rela- 
tively large errors. The results in the E* and S* chan- 
nels are reported for the first time. Some of our lattice 
result can be considered as predictions in cases where 
the spin-parity assignment is unknown in the PDG. The 
spin-l/2± states extracted from the spin-3/2 interpolat- 
ing fields are in good agreement with those from the stan- 
dard spin-1/2 interpolating fields, providing a non-trivial 
check of the calculation. Furthermore, we find that the 
projected spin-1/2 states are quite different in the full 
and the truncated interpolating fields. This means that 
as far as spin-3/2 states are concerned, the first term of 
the interpolating field is sufficient. However, if the masses 
of spin-1/2 states are computed from such interpolating 
fields via spin projection, the full interpolating field is 
required. 

As an independent check of the calculation, we car- 
ried out a parallel calculation of the usual baryon decu- 
plet (isospin-3/2 and spin-3/2) on the same lattice with 
300 configurations, using the same projection techniques. 
The pattern in the A states is consistent with the previ- 
ous calculation 25] and with experiment. This reinforces 
the efficacy of the methods used in separating the spins 
and parities. The results in the S* and S* channels are 
new and are used to shed light on the spin-parity of some 
states in the PDG. 

Having established the signals and the methods used 
to isolate the 3/2± and l/2± states, improvement can be 
made in a number of areas in future studies. First, higher 
statistics (probably on the order of 1000 configurations) 
are needed to beat down the errors in the weaker spin- 
projected states. Second, smaller pion masses are desired 
to perform a chiral extrapolation and make better contact 
with experiment. Third, both the lattice spacing and the 
box size should be varied to assess possible discretization 
effects. In the long run, the calculations should be done 
with dynamical gauge configurations in order to assess 
the effects of quenching in this sector. 
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isospin-1/2 and spin-3/2 interpolating field in Eq. @ is 
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TABLE 11: Mass ratios for the four states in the isospin-3/2 and spin-3/2 family (baryon decuplet) after parity and spin 
projections. The numbers in brackets are statistical errors in the last digits. The last column indicates the time window from 
which the results are extracted. 



m^/mp 



A(l/2+)/iV 
A{l/2-)/iV 
A(3/2+)/A^ 
A(3/2-)/Ar 



E*(l/2+)/iV 
E*(l/2-)/iV 
E*(3/2+)/iV 
E*(3/2-)/7V 



S*(l/2+)/iV 
H*(l/2-)/iV 
H*(3/2+)/iV 
5*(3/2-)/iV 



0.948(1) 0.935(1) 0.916(2) 0.889(2) 0.847(2) 0.817(3) 0.777(4) 0.723(4) 0.643(5) 10-15 



1.35(33) 
1.17(2 
1.03(1 
1.16(1 



1.07(9; 
1.02(2 
0.89(1 
1.03(1 



0.96(8 
0.88(2 
0.75(1 
0.90(1 



1.38(34) 
1.19(3) 
1.03(1) 
1.18(2) 



1.42(36) 
1.21(4) 
1.04(1) 
1.22(2) 



1.47(35) 
1.25(5) 
1.06(1) 
1.26(2) 



1.53(36) 
1.29(7) 
1.08(1) 
1.32(3) 



1.57(35) 
1.33(8) 
1.09(2) 
1.36(4) 



1.62(37) 
1.37(10) 
1.11(2) 
1.41(5) 



1.68(37) 
1.42(13) 
1.14(2) 
1.46(6) 



1.76(40) 
1.47(18) 
1.18(3) 
1.54(8) 



7-9 
10-12 
11-13 
11-13 



1.12(11) 
1.06(3) 
0.92(1) 
1.08(2) 



1.18(13) 
1.11(3) 
0.95(1) 
1.13(2) 



1.26(15) 
1.17(5) 
0.98(1) 
1.20(2) 



1.35(17) 
1.26(7) 
1.04(1) 
1.28(3) 



1.40(19) 
1.31(8) 
1.07(2) 
1.34(4) 



1.47(21) 
1.37(10) 
1.11(2) 
1.41(5) 



1.54(23) 
1.44(12) 
1.16(2) 
1.49(6) 



1.63(26) 
1.53(17) 
1.23(3) 
1.59(7) 



7-9 
10-12 
11-13 
11-13 



1.03(10) 
0.94(3) 
0.79(1) 
0.96(2) 



1.11(11) 
1.02(4) 
0.85(1) 
1.04(2) 



1.20(13) 
1.10(5) 
0.91(1) 
1.13(2) 



1.32(16) 
1.22(7) 
1.00(1) 
1.25(3) 



1.39(18) 
1.29(8) 
1.05(2) 
1.32(4) 



1.47(21) 
1.37(10) 
1.11(2) 
1.41(5) 



1.56(24) 
1.46(12) 
1.19(2) 
1.51(6) 



1.68(27) 
1.58(15) 
1.29(3) 
1.65(7) 
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